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THE CRITICAL SHEAR LOAD OP REOTAlTaULAR PLATES* 

By Edgar Seydel 



!• SCOPE OF PROBLEM 



In the stress analysis of plates in shear the deter- 
mination of the critical shear load, i.e., the load at 
which an originally flat plate "begins to buckle, is of 
primary significance, because when this load is exceeded 
the plate breaks or at least assumes a different stress 
attitude (if the plate is very thin) • The analytical meth- 
od developed by Time schenko** and Reissnor*** for the crit- 
ical load in shear has been extended to include the orthog- 
onal- ani sotropio (or tho tropic) plate .**** 

The premises for the solution of the problem are: 
The perfectly flat, rectangular, thin, homogeneous plate 
of constant thickness is supported at the edges without 
being restrained. The edge supports are rigid and the con- 
nection of the plate with the edge supports is such that 
the plate along the four boundary lines cannot bend out of 
its plane. Along the four edges the evenly divided, shear 
load t of constant intensity, is applied, which is in 
outside equilibrium and which, provided the shear is so 
small that the plate remains flat, impresses the same shear 
stress in the entire plane of the plate. (See fig. 1.) 

In view of the use of plywood plates and corrugated 
metal strip or othor similarly stiffened plates, we shall 



♦ " Ausbeul-Schubl ast rGchtertkigfir Platter.,'* Z.F.M,, ?cl: - 
ruary 14, 1933, pp.„78-.33. 

**Timoschenko : "Uber die Stabilitat versteifter Flatten, 
Zeitschrift "Der Bisenbau," Vol. 12, Kos, 5 and 6, 1921, 
pp. 147^163. 

***Bergmann and Reissner: "Uber die Knickung von recht- 
eckiffen Flatten bei Snhnhne^nspruchnng. " Z.?,M., Vol. 23, 
Wo. If^ 1932, pp. 6-12. 

♦***"Uber das Ausbeulen von recht eckigen, isotropen oder 
orthogonal-ani sotropen Flatten bei Schubbeanspruchung, " 
(To be published in Ingenieur-Archiv, 1933.) 
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treat not only tho isotropic "but also tho orthotropic plate 
in which the rigidity in two directions at right angles to 
each other is different. Heroin it is ass-amed that those 
two directions (which, for instance, in the laminated plate 
corrosTDond with the dirocti on ox ij. 0 fi'bors in the outer 
and inner plies) are parallel to the plate edges. 



II. FORMULAS FOR CALCU1ATIH& THE CRITICAL LOAD 
in SHEAR OF RECTAH&ULAR PLATES 



1, Isotropic Plate 

Let a he the long side of the rectangiilar plate, 
b, the short side, and 6 the thickness. Material prop- 
erties as far as they enter into the discussion are char- 
acterized by the elasticity modulus E and the transverse 
elongation figure (Poisson's constant) v. Then the plate 
stiffness is 



12 (1 - v^) 
and the critical shear load is: 

''^-^ = °" l^sY 

with Cq^ as a coefficient contingent upon aspect ratio 

p = h/a to he read from figure 2, For the infinitely long 
plate (p = 0), Ca = 13.2 (or more exact, 13,165); for the 
square plate (p = l) it is Ca - 23,0; for other aspect ra- 
tios p, c^ lies between these two limits, 

Z. Orthogonal-Ani so tropic Plate 

Assume a and b are the two sides of the rectangu- 
lar, or thogonal-ani so tropic plate and 6 the plate thick- 
ness; further, let the direction of side a be called x 
direction, and that of side b* y direction. (Seo fig. 1.) 

Whereas the determination of the elastic form changes 
in the isotropic plate can be effected with the data on E 
and it requires in the case of the orthotropic plate 

the data on E and u for the directions given by the or- 
thotropic anisotropy and in addition the data on shear mod- 
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ulus. Both E and v are admittedly dependent on each 
other, so that one of these four quantities is already 
defined "by the three others. 

If Ex is the elasticity modulus for length changes 
of the fibers in the x direction and the correspond- 

ing transverse elongation factor, and further, if Ey is 
the elasticity modulus for length changes of fihers m the 
y direction and Vy the corresponding transverse elonga- 
tion factor, we have:* 

Ex = '^X Ey 

Now let & "bo the shear modulus. 

Then the stiffness per unit length is: 

/_-v „ „ jflexural stiffness in x direction 

\^o)x ^2 \(i5y flexure ahout y axis) 

riPT^ = TT ^ j f loxural stiffness in y direction 
^EJ^y £y | ^-^^ flexure about the x axis) 

4(GJ)xy - 4G 2,-. torsional stiffness & -g-y 

Then we compute 

(EJ)x 



Dp = 



1 - Vx ^y 



(EJ). 



1 - ^x ^y 



= i (^x Da + ^y Di) + 2(Gj)xy 

1 D3 ((reciprocal criterion of the ortho- 

^ " T^^dT'S" plate) 

If |- = 1, we compute 



*Luf tf aiirtf orscliung. Vol. 8, No. 3, p. 78, or D.V.L. Year- 
took, 1930, p. 242. 
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Tlien the critical shear load is: 

4 



*lrr = Ca 



The Cg^ factor in this formula is dependent on Pa 
and l/-^ and may be tal-en from figure 3.** 

But if 7 > If we figure \rith 
and the critical shear load is: 



*kr = 



(t/2r 



The definition of 0-5, which occurs in this forraula, 
is explained in section IV of the report cited in footnote 
**** (first page)* 

Plates with evenly disposed reenf or cement s can, unless 
the spacing is unduly grer,t, "be approximated as homogeneous 
orthotropic plates and (SJ)^^, (EJ)y., and 4(GrJ)2j.y ^^^e then 
expressed hy mean values. The method of computing this 
stiffness, say, for corrugated metal sheet, has "been de- 
scribed in D.V.L, report xTo» 230, "Coliimn Tests in Shear 
with Corrugated Metal Sheet." (D.V.L. Yea,rbook, 1931, pp. 
233-245.) 



*If should yield a figure greater than 1, the nota- 

tions for the two rectangular sides should be exchanged, 
i«e», if the longer side was denoted by a, it should then 
be changed to b and vice versa. The same applies, of 
course, to and . Ordinarily the long side is called 

a, and the short side b, but it may happen that the stiff- 
ness in direction of the long side exceeds that in direc- 
tion of the short side, in which case then the switch rau.st 
b o mad e • 

**In this diagram we chose the reciprocal criterion 
in preference to ^ as parameter, to ensure linear inter- 
polation for any criteria (^'ithin range of 1 ^ ^ ») • 
The curves Ca (?a) were obtained by an approximation 
method. (See footnote (first page.) Although it is 

likely that the so ctilculated values do not substantial- 

ly deviate from the exact values. 
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III. EXAMPLES 
Isotropic Plate 

!• Example: 

Duralumin plate, 0,3 mm thick. 

B = 7.5 X 10^ kg/cm^ v ^ 0.25 

a) Dimensions: length 45 cm, width 15 cm 

B = O! = 7.5 X 10^ X 0,02^^ 0^533 j^^/^^ 

12 (1 - V^) 12 (1 - 0.25^) 

Por ^ = ~ = = 0.333, figure 2 gives: 
a 

Ca = 14,3, 

that is, a critical shear load of 

= 0.13& kg/cm 

b) For an identical hut square plate having sides 
a - "b == 15 cm (i.e., p = 1) , we have: 

Ca = 23.0 

and consequently, 

*kf = 23.0 X = 0.218 kg/cm 

7.5 

Orthotropic Plate 

2. Example; 

Plywood plate, 2 mm thick. The elasticity modulus* 

*The elasticity factors really should "be determined hy 
flexure and torsion tests. See Hertel's report on the 
shear moduli of laminated and ply wood, D.V.L. Yearbook, 
1932, III, pp. 43-52. 
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in the direction of the i'^rain (i.e., of that of the outer 
ply) is assumed at 150,000 kg/cm^, the corresponding 
transverse elongation factor at 0«4, the elasticity modu- 
lus across the grain (i«e,, of that of the inside ply) on- 
ly 15,000 kg/cm^ one-tenth of 150,000), the correspond- 
ing transverse elongation factor likowiso 0,04 (i»e«, one- 
tenth of 0,4), and the shear modulus G == 10,000 kg/cm^. 

a) Dimensions of plate: 40 hy 60 cm; direction ox 
grain parallel to short (40 cm) side, 

a = 60 cm Ex = 15,000 kg/cm^ = 0,04 

b = 40 " Ey = 150,000 « = 0.4 

& = 10,000 " 

Dx = ' ^2 J = 3^ J 15,250 J kg,W, 

1 - ^y 0,984 ' 

Ev 150000 , ^ 

^2 = r~ ~ ^ " ^^57984 ^ " 152,500 J kg/cm^ = 10 . 
y 

D3 = Uy Di + 2 G J = (5,100 + 2 X 10,000) J=26,100 J kg/cm^ 

on 

= 0.54, 



- P3 _ S610 0 



^AD^ '1)2 15250n,^0 



t */ Di 40 1 



Pa = - / — = — X = 0,375 

According to figure 3, c^ - 11*8 

t/^T^V^^ 15250 y 10^ ^3 

*kr = ^a " = 11-8 20^ X 0.367 X 10 

*kr " 1«69 kg/cm 



b) How assume that the plywood plate is divided into 
four equal parts "by three strips arunning i,n the direction 
of the side, '.vhich is 40 cm long (the ni.c.yt side). These 
stiffening strips (as well as the edges) shall "be s.o stiff 
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as to prevent "bending out of the plane of the plate when 
the latter huckles; the plate to ho supported without fix- 
ity at the stiffening strips as well as at the edges. The 
grain to run parallel to the long (40 cm) side* 

7" = 0t54 (unchanged) 
*y — r- 

Pa = i§ = 0-667 

According to figure 3: 

ca 13.5 



15250 ^/lO ^-3 
tw = 13 — 2 ~ X 0.567 X 10 

*kr " 4.35 kg/cm 

The stiffening strips raise the critical shear load 
approximately 2.6 times, 

c) The ,same example as a) , hut with grain running 
parallel to the 60 cm side. 

Ex and Ey (Di and D2 ) would then have to be switched, 
but then p^^ > 1, so we put 

According to figure 3: c^ = 15.4 



_ . 15250 / 10^ V n ^^r, V in-^ 
iv^.u 30^" " — ' 

tkr = kg/cm 



In this case the critical shear load amounts to 0.58 
of the value of example a) • 

d) The same example as b) , but grain as in example c) 
parallel to 15 cm (respectively, 60 cm) side. 
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Ca = 11.3 

t 11,3 liaSS^pZ X 0.667 X 10-= 

7.5^ 

*kr ~ kg/cm 

With the grain running in this direction, the stiffen- 
ing strips augment the critical shear load appr oximatels'' 
11,7 times, "by changing the direction of the .^rain against 
example h) to about 2.5 times, 

3, Example: 

Sqiiare, corrugated duralumin strip: a = "b = 90 cm 
Elasticity modulus of strip: E = 750,000 leg/cm^ 
Transverse elongation: v = 0.25 

Cross section of corrtigated sheet: 21 = 3,0 cm 

(•length of 
corrugation) 

2f = 0.915 cm 

(height of 
corrugati on) 

Thickness: 6 = 0.039 cm 

Applying the formulas given in report ITo. 230 of the 
D, V.l/ learhook, 1931, pp. 233-245 (table II and figures 
5 and 6) , we find: 

Di = 3,27 kg/cm 

= 3,790 kg/cm 

D3 = 3.585 kg/cm 
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= 8.3 



=8.3 




= 2.7 kg/cm 



Although the plate is square the coefficient Cg^ is 
only about 2 percent higher than Ca = 8.125 resulting 
for an infinitely long orthotropic plate (Pa = 0) ^itlt 



At the recommendation of and in collaboration with 
Professor Reissner, we made a series of shear tests on 
thin duralumin sheets (occasionally also with plates of 
German silver, tin and celluloid), in which the theoreti- 
cal assumptions were to he fully corrohorat ed« In order 
to ensure a free edge support in the experiment the hori- 
zontally rigged up plate was restrained on top and bottom 
side by knife-like edge rails. These rails formed a rec- 
tangle, one side of which was 15 cm long, whereas the 
other side was 1, 2, and 3 times as long. The test plate 
itself was slightly greater than the rectangle formed by 
the rails, extending about 1.5 cm beyond each side. About 
0,5 cm distant and parallel to the rails we rigged up four 
loading rails hinged by bolts and provided with pins, 
spaced about 1.4 cm, which fitted the corresponding holes 
on the test plate. The loads which were mutually in equi- 
librium were applied at two diagonally opposite bolts of 
the hinged quadrangle, through which the loading evenly 
distributed, was initiated in the test plate. In order 
to disturb as little as possible the irregular elongations 
normal to the edge which result from buckling, the holes 
in the plates were elongated. 

The buckling was recorded by a Zeiss dial gage. Since 



NOTE.- The metric values used in this T.H. can be convert- 
ed to the English units by use of the following 
conversion factors. 



IV. 



EXPERIMEHTAI DATA 




cm3 X .061023 = cu.in. 

kg X 2.20462 = lb. 

kg/cm X 5.59977 = lb. /in. 

kg/cm2 X 14.2235 = Ib./sq.in. 
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the pressure exerted "by the spring inside of the clockwork 
would signify a minor tut still appreciable bending load 
on the test plate, the spring tension was so balanced that 
the tip of the dial gage, hence the pointer in every posi- 
tion which the tip assumed, stood still owing to the in- 
ternal friction of the dial gage, but upon slightly tap- 
ping, moved a little in the sense of the spring force ef- 
fect. Prior to each reading the tip was first set slight- 
ly lower than the test point (metal plate) » then the tap- 
ping was continued until the tip touched the test point. 
This contact was recorded by a circuit and a galvanometer 
deflection* 

To illustrate the test results the deflections of the 
plate at a test point (a point at which a maximum buck- 
ling was anticipated, generally in the plate center) were 
plotted against the progressively Increased loadings. In 
addition we plotted contour-line charts based upon the read- 
ings of a greater number of test stations for each l«ad 
stage, which brought out the differences in the buckling 
at different load stages (generally the differences between 
a load beyond the theoretical buckling load and zero or 
very small load). 

Theoretically the aspect of buckling at a test point 
when shown versus the load should be as the dot-dash line 
in figure 4. There is no sign of buckling as the load in- 
creases till the critical load is reached. Then buckling 
occurs suddenly, i.e., the curve should at its start (on 
the abscissa axis) have a tangent perpendicular to the ab- 
scissa axis. The experiments failed to reveal such a 
curve as shown in figures 4 to 7. Even the measured con- 
tour maps (figs. 8 to 10) show a not unsubstantial dis- 
crepancy from those obtained according to the theory.* 
Said discrepancies in theory and test reveal that the the- 
oretical aqfsumptions were not fully complied with in the 
experimental. In particular, in the tests the effect of 
even very ^mall initial buckles appears perceptible, which 
could not be altogether avoided. 



*The theoretical buckling areas for plates of 1:1, 1:2, 
and 1:3 aspect ratio are given in report cited in foot- 
note *** (first page) • Subsequent but more accurate cal- 
culations revealed that the joints (w = 0) of the buck- 
ling surfaces should Lave, where they go over into the 
edge lines, tangents at right angles to them. 
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We likewise experimented witli corrugated metal plates 
whicht ^ith certain approximation, may be considered as 
orthotropic plates. Several of these tests are described 
in the report cited, No. 230, D.V.L, Yearbook, 1931, pp% 
233-245; other similar test s, not touched upon in this re- 
port, yielded substantially the same results. In these 
experiments made in connection with the theoretical treat- 
ment of infinitel y lo ng orthotropic plate strip the value 
for Pa ^ b/a\/'l3^75^"^ lies between 0.034 and 0.072t The 
validity, surmised at that time on the basis of simple 
considerations, to introduce the value Cq, for such small 
Pa as afforded by ^q, - 0 (infinitely long plate) is now 
proved by the c^ (^q) curve in figure 3, at least for 
the freely supported plate. Admittedly the plate edges 
in these tests were (almost rigidly) restrained (divergent 
from the problem treated in the present report). 

The set-up for these experiments with one or two 
square corrugated duralumin sheets was somewhat differ- 
ent. The sides of the sheets were 90 cm long* The sheets 
consisted of three strips riveted together and had (apart 
from slight deviations definable only by the most accurate 
measurement) the same, sectional dimensions. Four duralu- 
min tubes formed the four edges, to which, on each side, 
a sheet was riveted. In order to initiate the shear more 
evenly through the individual rivets into the sheets, the 
resultant load allocated to each edge support was applied 
in the center (instead of at the end) , 

The riveting of the corrugated sheet produces an elas- 
tic restraint on the edge supports, which, compared to the 
free support, effects an upward shift in critical load. 
The buckling at different points was measured in the me- 
dian line of the sheet perpendicular to the direction of 
the wave ridge. Whereas the majority of points revealed 
from the start certain buckling, at first gradual, then in- 
creasing, the course of buckling at one point, versus the 
load, corresponded throughout to the course anticipated 
by the theory, as shown in figure 11 • The critical load 
must lie between the figures for free edge support and 
those for rigid restraint, (The value for rigid restraint 
in the figure corresponds to that for the infinitely long 
plate, but might be only slightly too small for the test 
sheet, because p^ 0^171,) 

The first test actually revealed no appreciable buck- 
ling up to a certain load (critical shear load). This 
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critical load lies in the zone imposed by the elastic 
edge restraint. Only minor permanent form changes re- 
mained after unloading, which, however, effected certain 
deviations in the course of the curve when the experiment 
was repeated (as likewise anticipated, according to theo- 
retical considerations). When repeated after having been 
previously uiiloaded, this deviation became still greater. 
The sheet failed at 2.3 times the buckling load, which 
was determined from the course of this curve. 



V. SUMMASI AND REPEREITCE TO 
THE PRACTICAL ASPECT OF THE PROBLEM 



The report adduces formulas for analyzing the critical 
shear load of a. simply supported square, isotropic, or or- 
thogonal anisotropic plate. Although arrived at by approx- 
imation method, these formulas should nevertheless be am- 
ply accurate for practical purposes, especially since the 
premises of the theory (unrestrained support at the edges, 
perfectly flat plate) never aire precisely fulfilled. The 
experimental data reveal that in very thin plates the dis- 
turbing effect of incipient buckling is quite considerable. 

If the plates are very long, rigid edge restraint can 
also be taken into consideration, by means of the formulas 
given in the D,V.L. Yearbook for 1930, p, 237, 

The interpretation of the strength of a plate in shear 
according to the ratio of plate thickness to plate width 
with a view to the problem in point here can be effected 
for three cases: 

1) That the plate is so thick that failure occurs 
by exceeding the shear strength, while the plate remains 
flat without signs of buckling; 

2) That the (pure shear) stresses in a thin plate 

up to attainment of critical shear load are relatively low. 
Owing to the buckling, however, due to this load, rapidly 
increasing flexural stresses are induced which may after 
exceeding the critical load lead to failure of the plates 
in which case, the critical load is of primary signifi- 
cance for the strength. 

3) That the flexural stresses contemporary with the 
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buckling are comparat iveljr low in very thin plates, such 
as used in aircraft construction; that they remain (apart 
from harmless purely local esicesses) "below the yield lim- 
it, so that, finally, the "break occurs "by exceeding the 
tensile strength after forming the Wagner tension diago- 
nal field. Starting with the critical load the elastic 
form changes occurring in the sense of the applied loads, 
increase in greater measure than "before. Otherwise, the 
critical load has no significance for the evaluation of 
the strength. 

Admittedly, no sharp line can "be drawn "between cases 
2 and 3. Por the "bar stressed in compression, case 2 is 
identical with the Euler load (column failure), case 1, 
with a short thick column which fails before reaching the 
stability limit by exceeding the compressive strength of 
the material, and case 3, with a very slender bar which 
can be loaded until both ends meet and even extend beyond 
without inducing permanent form changes, although this 
case is of scarcely any practical import for the compressed 
bar, in contrast to the plate stressed in shear. 



Translation by J. Vanier, 
National Advisory Committee 
for Aeronautics. 
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Figs. 1,2,3 

Figure 1. - Rectangular plate 

load in shear t, 
the grain direction running- 
parallel to the plate edges 
(as, say, of a plywood plate) 
.was used as "basis in the 
treatniant of the orthogonal- 
anisotropic plate. For load- 
ing "by Liomant Mi the flexural 
stiffness Di (elasticity mod- 
ulus E-) is decisive, for 
loading oy moment Ms the flex- 
ural stiffness D2(elasticity 
modulus K^) is significant. 



Figure 2. - Coefficient 

tkr (W2)2 



of the critical shear load 
of an isotropic plate versus 
aspect ratio p=lD/a 




Figure 3. - Coefficient 

ca=:feJM2if 

\/%i>2^ 
of the critical shear 
load of an orthogonal - 
anisotropic plate versus 

_ D ^/^i .... 



versus 
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Figs* 4.5,6,7 




flZTt 



3 6 9 12 15 kg 
P 

Figure 4. - Test data: "buckling w 

(in 10***''^ nm) versus 
load P ( kg), recoi-ded in center of 
0,2 mm gage square dural plate 
(a= t= 15 cm) by loading and un- 
loading. ? and Pi^^f (theoretical 
critical load) are the resultant 
loads falling into tlcK-: plate diag- 
onals. The dot- daslied line shows 
the anticipated curve under ideal 
test conditions. 
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Figure 5.- Test data: eame 
as in Fi^re 4, 
except that "b/a = 1/ S (13 = 
15 C5m, a = 30 cm) 
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Pig-are 4, except tliat 
■b/A= 1/3 ( 1)= 15 cm, a=45cm) 




0 30 60 SO 120 150 lig 
P 

as in Figure 4, 
except plate thicleiess, 
0.428 irni, b/a= I/3 ( h = 15 
ciu, a= 45 cm) 




Figs. 8,9,10 



Figure 8.« Test data; contour lines of Truckling surface of a C.3 
™i gage Square dural plate at P= 50 kg (= 3.2 Pkrf ) 




Figure 9.- Test data! contour map of a 0.2 mv. gage dural plate 

(Va= 1/2) "by P=: 42 kg (= 8.47 Pj^f ). (Discrepancy of 
buckling versus P = 2 kg) 




\< a=: 3lD= 45 cm ^ 



Figure 10,- Test data: contour map of a 0.2 vm gage dural plate 
(Wa= l/3)ty P= 30 kg. (=4.67Pwf) 
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Fig. 11 



mm 




0 100 200 300 400 500 600 700 kg 
P 

Figure 11.- Test data; tuckling w versus load P (resultant load 

in diagonal direction) , recorded on a square cor- 
rugated sheet ( side = 90 cm section: corragatian heigth / cor- 
rugation length / sheet thickness ( in mrn) = 9.15/ 30, O/ 0.39) 
"by maxltiple loading, (i^crf and rj^re ) i<^iiote the theoretical 
critical load hy simple support and rigid plate edge restraint. 
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